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Abstract 

Recently it was suggested that the observation of superluminal neutrinos by the 
OPERA collaboration may be due to group velocity effects resulting from close-to- 
maximal oscillation between neutrino mass eigenstates, in analogy to known effects in 
optics. We show that superluminal propagation does occur through this effect for a 
series of very narrow energy ranges, but this phenomenon cannot explain the OPERA 
measurement. 



1 Introduction 



Recently the OPERA collaboration reported a measurement of the average time taken for 
neutrinos (un up to % level contamination) created at CERN (CN) to arrive at the Gran 
Sasso Laboratory (GS) compared to the time taken travelling at the speed of light in vacuo 
(c). They found an early arrival time of approximately St = 60 ns, which corresponds, at a 
significance of 5.0cr, to faster-than-light travel by a positive fraction Sv = v — 1 ~ 2.37 x 10~ 5 
PQ. (In this paper we set Planck's constant and the speed of light in vacuo h = c = 1.) 

The OPERA result has inspired many papers (for a selection see [21 [3l 01 [5j [6] ) . Here we 
focus on the one by authors Mecozzi and Bellini [3]. They have suggested an interesting pos- 
sible interpretation in one-to-one analogy with the established dispersion and interference 
effects found in Mach-Zehnder interferometers and optical fibres with polarisation mode 
dispersion, where a superluminal group velocity has been both predicted and measured 

m-m- 

This effect does not imply that a signal is exchanged at faster than the speed of light, 
in violation with special relativity. Instead it arises from constructive and destructive 
interference deforming the leading and trailing edges of the pulse. 

Mecozzi and Bellini do not analyse whether such superluminal propagation is possible 
in practice in OPERA, and indeed in neutrino experiments in general. In this paper we fill 
that gap. We show that remarkably such an effect does occur in the OPERA experiment, 
causing superluminal 'spikes' in the neutrino velocity at various 'critical' values of the 
neutrino energy inside the neutrino beam energy spectrum. However this behaviour results 
in the wrong energy dependence and, once averaged over the full energy spectrum, is far 
too small to explain the OPERA result. 

Even though the spikes of superluminal group velocity caused by near-maximal neutrino 
oscillations cannot explain the OPERA measurements, it is an effect that is nevertheless 
interesting in its own right. We analyse carefully the size of the effect and the conditions 
that are required to realise it. In particular, we point out that the value of the maximum 
achievable superluminal velocity depends crucially on corrections which were neglected in 
ref. [3]. We will see that these corrections recover the result one would expect intuitively: 
neutrinos do not escape the domains of their mass-eigenstate wave packets, which travel at 
the expected subluminal speeds, therefore the maximum increase in speed is entirely due 
to the quantum mechanical uncertainty in position divided by the time of flight. 
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We note that this effect has by now been treated in two other papers. In ref [3], it was 
pointed out that the OPERA experiment might be understood as a manifestation of weak 
measurement. This is in fact another way of interpreting the same effect however, as in ref. 
[3], there is no attempt to determine whether this is realised for neutrinos in practice. 

Independently, the weak measurement interpretation was also treated in a paper that 
appeared after ours [5], where also a numerical estimate is provided. This is in agreement 
with our more detailed calculations, however the limiting corrections, multiple peaks and 
energy spectrum are not treated. On the other hand the treatment of the effect is performed 
in position space with Gaussian wave packets, which exposes particularly clearly the root 
cause as being through the small relative displacement of the mass-eigenstate wave packets 
and their near-cancellation. 

We will see that in order to maximise the effect we will need a mixing angle such that 
sin 2 (29) is as close to 1 as possible. It is known that 6±2 ~ 0.86, however the effective value 
of sin 2 (2^12) in rock is close to zero [ID]. #13 is tiny [IT], so we are therefore left with #23 
for which only the limits 0.92 < sin 2 (2^23) < 1 are so far known [10J. #23 drives O v T 
mixing and takes place practically as in vacuum [TDJ. Clearly then we can work effectively 
with 2- neutrino mixing between z/„ and v T . (Actually, this is not quite true: we know that 
the other mixing matrix parameters will supply corrections in exceptional circumstances. 
This will be discussed in the conclusions.) 

The rest of the paper is organised as follows. In the next section we review the derivation 
in ref. [3], including the notation and adaptations we will need to match the OPERA 
experiment. In particular we take care to include corrections neglected in ref. [3] but which 
are crucial in limiting the maximum size of the effect. To carry through the derivation, 
we will need to make some assumptions about the coherence of individual neutrino wave 
packets. In section 3 we show these are correct. In section 4 we show that while there is an 
effect, it is far too small and has the wrong energy dependence to agree with the data. 

In section 5, we draw our conclusions. 

2 Review and development of ref. [3] 

Our treatment follows closely ref. [3] but with adaptations and corrections. We write the 
mass eigenstates in standard convention [10] as: 





(1) 
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\p,u 3 ) = sin(0 2 3)b,^) + cos(6 23 )\p,v T ), (2) 

where p is the momentum of the neutrino in the beam. For momenta much larger than the 
masses, which is the case here (since the masses are known from tritium decay to be < 2 
eV [TO]), energies are given by 

E 2 = E + E/2 and E 3 = E - E/2 (3) 

where 

Eo(F) = p+ a^4^ p+ m and m = = (4) 

We start at CERN with a normalised pure ia, wave function |^o) a t time t = 0, and 
position xo = (V'oM^o)- Projecting on mass eigenstates using ^ fi){p', we have at 
later times 

|V> t ) = fdk(k\i; )e- lEot { S m(e 23 )e' tEt / 2 \k,u 3 ) + cos(0 23 ) e lEt / 2 \k,u 2 )}, (5) 



where ipo(k) = {k\ipo) describes the spread of momenta in the initial neutrino wave packet. 

To find the wave function when measured at Gran Sasso, we collapse \ipt) using 
(p 1 , Vfj\ to obtain 

4>'(p') = Mp'W) (6) 

where 

F(p') = { e- lEt ' 2 sin 2 9 23 + e iEt ' 2 cos 2 6 23 ] e"^ 4 , (7) 

and we have replaced the initial ipo(p') with an effective measured ipi(p') to allow for further 
decoherence in the momentum space of the particle during its measurement. 

We need to normalise ([6]), which requires dividing by fdp'\ip'(p')\ 2 . We now assume, 
as effectively was done in ref. [3], that ipi(p') is strongly peaked around the momentum 
p' = p, and that the variation in |.F(p')| is gradual in comparison. We will confirm these 
assumptions in the next section. We note that in this case Jdp'\ip'(p')\ 2 w \F{p)\ 2 and so 
the Gran Sasso wave function is typically normalised to good approximation by dividing 
by |.F(p)|. This approximation breaks down if p takes a value such that F{p) is close to 
vanishing (as is clear because fdp'\ip' (p')\ 2 is positive definite since F cannot vanish for all 
momenta). We will return to this in the next section. 

For now we take the expected position of the Vn at Gran Sasso, to thus be given by 

{x) ^\FW 2 l dP ' r ^ )i ^^' {P% (8) 



which depends on the time of arrival through F. Substituting (|6|), we have that the integral 
is given by N\ + N 2 where 

iVi = Jdp'\F\ 2 iP* 1 i-^ip 1 and N 2 = j dp' \^\ 2 F* i-^ F (9) 

Using again the fact that ipi is strongly peaked we see that 

N 2 ^F*(p)i^-F(p). (10) 

Since (JSJ) is the expectation of an Hermitian operator it must give a real answer. There- 
fore the imaginary part of the above must get cancelled. This can be seen to be true by 
integrating Ni by parts to get: 

iVT = + Jdp'\^\ 2 i^\F\ 2 . (11) 

and thus extract the imaginary part. The remaining real part of N% is approximately 

\F\ 2 (p) Jdp>rii^4>i~\F\ 2 (p)x . (12) 

Putting it all together we see that 

(x) = x + v g t « x + ^yp^ft |F*(p) i^- F(p)\ , (13) 

where we have identified the second term as t times the group velocity. Finally, substituting 
([7)1 and using we have 

\F(p)\ 2 = 1 - sin 2 20 23 sm 2 (\Am 2 23 \t/(4p)) (14) 

and evaluating the numerator: 

v g = 1 + 5v d + 6v sup (15) 

where 



*>* = -^r do) 



(^23 ) 

2p 2 

is the expected negative 'classical' correction to the speed of light as a consequence of the 
average squared mass of the two mass eigenstates, and the last term is a correction resulting 
from interference between the two mass eigenstates: 



Am| 3 cos 2#23 

" 4p 2 1 - sin 2 20 2 3 sin 2 ( | Am| 1 1 / (4p) ) ' 



Sv « 23 fgg fi7l 

UUsup~ . n __2n/l -2/1 A 2 \\' V X V 



This is the main result of ref. [3] (up to correction of the power of sin 2023 an d the explicit 
statement of approximation). 



3 Limitations to coherence and the strongly-peaked wave 
packet approximation 



Let us revisit the assumption that the momentum of the measured neutrino wave packet 
ipi(p') is strongly peaked about p' ~ p. This was used to derive the key formulae in sec. 2 
below ([7]), by assuming that the variation in |_F(p')| was by comparison more gradual. It 
should be clear that here we are not discussing the momentum spectrum of the ensemble of 
neutrinos in the beam, which is very broad — as we will describe in the next section, but 
rather the inherent quantum mechanical uncertainty in the momentum of a given neutrino 
as it is measured at Gran Sasso. 

We start with the neutrino wave function tpo(p') produced in the CNGS decay tunnel 
at CERN. The momentum uncertainty can certainly can be no smaller than that set by 
At ~ 5 ns, the smallest time features in the proton bunch [lJ^This corresponds to cAp = 
l/At « 1.3 x 1(T 7 eV. 

However, even if the proton beam is coherent at this level, it suffers decoherence on its 
way to becoming the z/« beam. Firstly, the protons impact the graphite target, producing 
the mesons (mostly pions) that will decay to muon neutrinos. Initially these mesons are in 
a quantum state together with the other products of the collision (including various nuclei) , 
however they then suffer decoherence from thermalisation in a hot target, both directly and 
also through their quantum mechanical coupling to the decay products. Assuming a target 
temperature of, say 300°cJ^| this limits the energy-momentum resolution to ksT ~ 0.05 eV. 

Finally the mesons decay in the decay tunnel and here further decoherence takes place, 
again through coupling to the decay products (in this case the muon) . Consider for example 
the decay of a ir + to /i + z^. The resulting quantum state takes the form: 



where Vv(q) is the wave function of the erstwhile pion, q and p are 3-momenta, and M 

stands for the matrix element for the decay. The muon is absorbed by a combination of rock, 

a Hadron stop and two muon detectors [U [12]. This allows the experimenters to measure 

the transverse coordinates of the proton beam spot when it hits the target to a precision 

1 In the new version of the OPERA experiment At ~ 3 ns, as set by the width of the bunch. 

2 This temperature is confirmed as the average maximum temperature in the revised version of the OPERA 

experiment [T|. 




phase space 




(18) 
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of ~ 50 - 90 /im. [12], however it is reasonable to assume that the rock itself localises the 
muons at the //m level (similar to emulsion — see below), even if this is not recorded. This 
corresponds to a momentum decoherence of order 0.2 eV/c. Through (|18p this decoherence 
is transferred to the neutrino. 

We conclude that the chief limiting factor on the coherence of the initial neutrino wave 
packet is through the decay of the mesons and results in a wave packet with momentum 
uncertainty w ~ ±0.1 eV/c. 

For muon neutrinos that interact in the Gran Sasso detector, the impact spot is dis- 
cernible in the emulsion at the ^m level [13] ; we can expect similar localisation in the rock 
in front of the detector for the external events. Therefore the act of measurement results 
in a momentum spread in the wave packet of similar size to that in the initial packet. 

We now contrast this with the variation in |_F(p')| from ([7]). We see that this is controlled 
by the phase 

m = m = J*^, (id) 

using (|3J), where t can be taken to be L/c and 

L = 730,085 m (20) 

is the corrected distance from the average meson decay point in the CNGS decay tunnel to 
the origin of the OPERA detector pQ. At the average energy p = 17 GeV/c of the neutrino 
beam [1], and using Am^ ~ 2.43 x 10~ 3 eV 2 [10] . we have that has magnitude 0.26, and 
thus we see that |_F(p')| is indeed slowly varying even on changes of order Ap' ~ GeV. 
More generally we require 

2 



d(f> 




dp' 





30.0 keV\ 

< 1. (21) 



pi 

The minimum energy-momentum p' we can consider is set by the resolution of the Gran 
Sasso detector, which is E m i n « 1 GeV [13] . or more optimistically by the lowest threshold 
energy for the charge current interactions used to detect the v^. This occurs for the process 
v^n — > [v~p, which would imply E m { n > + ^m 2 t /m n = 112 MeV. Either way we see that 
the condition ([21 j) is amply satisfied. 

Although this establishes that is slowly varying compared to the fundamental 

uncertainty w set by the width of the wave packet, as we mentioned in sec. 2 we cannot 
trust (j!7p if F(p) is close to vanishing. This leads to extra conditions that must be satisfied 
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if we are to trust the approximation. From ()14j) . we see that takes its minimum 

value, cos 2 2023, when 4>{p) = ir + 2na, where a = 0, 1, 2, • • •, This corresponds to 'critical' 
values of momentum 

. . lAm^lt 1 1.43 ^ rr; 

P = *H,(.) = S^J^ = GeV/c. (22) 

Taylor expanding about these values we have 



j F(p')\ 2 = cos 2 2023 + -^-(1 + 2a) 2 sin 2 20 23 [Ap'] 2 + 0([Ap'f), ( 2:5) 



where Ap' = p' — p cr it- Thus from ([6]), at p = p cr u, we can compute the first correction 
to jdp'\ip'(p')\ 2 « |F(p)| 2 , from the finite size of the wave packet. It is independent of the 
details other than through the RMS width w 2 = ( , (/'i|[Ap / ] 2 |'0i): 



2 2 

dp'W{p')\ 2 = cos 2 2023 + ^^(1 + 2a) 2 sin 2 20 23 + 0(w 4 /pt rit ), (24) 



■ r 2 

Thus if 



7T It) 

|cot20 23 |»- (l + 2o) (25) 

^ Pcrit 

the first term in (|24p dominates and J'aVlV^p')! 2 ~ I-^(p)| 2 * s always a good approximation. 
If (|25p is not satisfied then p cannot be taken too close to p = p cr it(a): we can only trust 
the expression (fT7|) for 5v sup providing 

2 2 

|F(p)| 2 > cos 2 20 23 + ^--^-(1 + 2a) 2 sin 2 20 23 . (26) 

4 Confronting the experiment 

Given the average 17 GeV energy of the beam, the classical correction to the speed of light 
([To]) is negligible and can be neglected. For most values of p, Sv sup is also clearly negligible. 

However, (|17p provides a correction that is maximal at the critical values of momenta 
(|22p . Note that the requirement of detection at Gran Sasso sets the upper limit to a as 

a < ' 231 - -. (27) 

47T±!/ m j n Z 

This implies a can only be zero if E m i n = 1 GeV, or a < 6 for E m i n = 112 MeV. 

If 023 is infinitesimally close to the maximal mixing value of vr/4, then (|17p actually 
diverges at the points p = p C ritip). To see this more clearly, let 

2023 = g =F 56 and P = Vcrit + e, (28) 



where the sign is opposite to that of Am^, 59 > is a small shift from maximal mixing, 
and e is a small deviation in energy from the critical value. Then (|17p becomes 

|Am| 3 | 59 

sup "WwTA«j' (29) 

However, at e = 0, from (|25p . this is only trustworthy if 

\50\ > ~ — (1 + 2a) = 1.10 (1 + 2a) 2 x 1(T 10 (30) 

Pcrit 



Substituting this in ()29p and using ()22p . we see that the maximum super luminal effect that 
can be predicted is 

5v sup < — = 2.70 x 10~ 12 . (31) 
wt 

This formula has a simple intuitive explanation. 1/w is the spatial width of the neutrino 
wave packet. Therefore the upper limit on 5v sup is nothing but the maximum increase in 
speed that could result from the quantum mechanical uncertainty in the position of the 
neutrino. We see that in order to explain the OPERA measurement we would have to 
assume that the width of the neutrino wave packet created at CERN is 1/w 3> St = 60 ns 
which is not compatible with the analysis in the previous section. 

We note also that OPERA repeated the analysis concentrating on only those charged 
current events occurring in the OPERA target (where reliable energies could be measured). 
They split this sample into bins of nearly equal statistics, taking events with energy higher 
or lower than E me d = 20 GeV. With a significance of greater than 3a, they still see a 
superluminal velocity in the higher energy sample of the same magnitude (within errors) 
PP. If this group velocity effect was the explanation then we ought to find that the effect 
goes away in the higher energy bin. 

Finally let us show that, taking into account the spread in energies in the neutrino 
beam, the superluminal contribution from this effect is in fact some 12 orders of magnitude 
smaller than (|31|) . Only a very narrow spread in energies in the neutrino beam contributes 
significantly to (|29p . We see that the energies contributing must lie in the range 

p = p C rit ± e, with e < 2p crit 59/-K (32) 

In fact, from the talk given at CERN [T3] and earlier OPERA analysis [15], we can see 
that the energy spectrum is broad, rising approximately linearly from E m i n to a maximum at 
E = E me d, and then falling with a large tail reaching energies of ~ 200 GeV. Let the fraction 
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of neutrino events measured between momenta p and p + dp be n(p). If we model n{p) = pp 
for p < E me d, where p is a constant, then, since about half the events are found below E mec [, 
we have by p « l/E 2 med . It follows that n(p crit ) « Pcrit/ E 2 med = 3.6 (1 + 2a)" 1 x 1(T 12 eV" 1 . 
Since (fTT|) is sharply peaked around p = p C Hti the average contribution from ([29]) is to very 
good approximation 

fdpn(p)v sup (p) = ~|AmI 3 |£^U ( 33 ) 

where A is the number of peaks (1 or 7 depending on E m i n ). This is numerically 3^4 x 10 -24 , 
which is of course too small to measure. It is also competitive with the classical term (|16D . 
Using the same model for n{p) we have for the classical term 

fdpn{p)v cl {p) < ~^In , (34) 

where the inequality results from neglecting the contribution of energies greater than E me &. 
This term will dominate unless there is a large hierarchy between the two masses, in which 
case this is of the same order, making the detected average group velocity slower (faster) 
than the speed of light, depending on whether the larger (smaller) value of E m i n is used. 



5 Conclusions 

We have seen that if the mixing angle #23 is close to maximal then thanks to intereference 
between the two mass eigenstates, the group velocity for the muon neutrino wave packet, 
Q15p using P^|) and (fTT|) . develops a series of very sharp peaks at critical momenta 

\Am\n\t 1 ,„_v 

P = ^riM)= ^ — , (35) 

where a is a non- negative integer. It is surely interesting in its own right that super luminal 
propagation of neutrinos is possible in principle through such an effect. In the OPERA 
experiment this corresponds to energies 1.43/(1 + 2a) GeV and in practice a is bounded 
above by the fact that for sufficiently large a it is no longer possible to detect the neutrino, 
definitely a < 6. 

However from (|3ip we have seen that the maximum extra displacement caused by this 
is limited above by Ax = 1/w, where Ax is the width of the wave packet and w is the 
width in momentum space. This translates into a maximum increase in group velocity of 
Ax/t, where t is the naively expected time of flight. Therefore we see that the effect is not 
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in conflict with relativity and causality, being limited to no more than that caused by the 
inherent quantum mechanical uncertainty in position. We have argued that the wave packet 
has a width of order 1 fim, as a result of decoherence from the detection (whether recorded 
or not) of the associated muons. Even if this is not the case we noted that thermalisation 
in a hot target would limit the coherence length of the pion wave packet, and thus also the 
neutrino wave packet to ~ 4/xm. This corresponds to a maximum superluminal correction 
of 6v sup ~ 10 -11 c and is therefore far too small to correspond to the effect seen by OPERA. 

We note that this also provides potentially a severe constraint on proposals to explain the 
OPERA measurement via superluminal propagation of only certain mass eigenstates (e.g. 
a sterile neutrino in extra dimensions [B]) since the OPERA measurement corresponds to an 
extra displacement by approximately 20 m, whereas a relative shift between mass eigenstate 
wave packets of more than a few fim is enough to destroy the coherence of the observed 
neutrino oscillations. 

We do not reach the upper bound, Ax/t, on the maximum extra group velocity unless 
the mixing angle #23 is closer to maximal than a term of order w/p, cf. eqn. ()30p . Such 
finite-width corrections are negligible everywhere except in this case, and then only when 
close to the group velocity peaks, where they matter because of they correct the near 
vanishing of the denominator term in (|17j) . There are some corrections to the numerator 
that become important similarly since it is also close to vanishing, although we did not need 
them for our analysis. In practice the neglected small corrections due to the other neutrino 
mixing parameters (5, #13 and the effective value of #12 in rock) could also have an effect 
when close to these peaks for the same reason. The exact behaviour close to the peaks is 
therefore the result of a competition between these corrections, and will depend on which 
one dominates. 

None of this matters for the measurements performed in the OPERA experiment how- 
ever since they measure an effect that is averaged over the broad energy spectrum of the 
neutrinos. The average effect would give a strong dependence on energy, which is not seen 
in the OPERA experiment, if it weren't for the fact that the result is anyway too small to 
measure. Furthermore we have seen that, once averaged over the energy spectrum of the 
neutrino beam, the resulting net correction (|33p is ~ 10" 23 , at best of the same order as 
the negative (i.e. subluminal) classical correction. 

It seems very difficult to imagine a scenario where this effect could explain the OPERA 
measurement. We note that the neutrinos with superluminal group velocity would be found 
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preferentially at the leading edge of the neutrino bunch. But from the plots [HE], it is clear 
that there is no tail of early arrivers: the shift due to an excess of super luminal neutrinos at 
the leading edge cannot be much more than the measured 5v, and besides we have too few 
of them by this effect in the energy spectrum. Similar comments apply to the superluminal- 
depleted population at the trailing edge. With the new measurements using shorter pulses 
[1] all such types of explanation are ruled out. 



Acknowledgments 

The author thanks the STFC for financial support. 



References 

[1] T. Adam et al. [ OPERA Collaboration ], "Measurement of the neutrino velocity with 
the OPERA detector in the CNGS beam," 
|arXiv:1109.4897 v2 [hep-ex]]. 

[2] D. Fargion and D. D'Armiento, "Inconsistence of super-luminal Opera neutrino speed 
with SN1987A neutrinos burst and with flavor neutrino mixing," arXiv:1109.5368 
[astro-ph.HE]; G. F. Giudice, S. Sibiryakov, A. Strumia, "Interpreting OPERA re- 
sults on superluminal neutrino," [arXiv: 1 109 . 5682 [hep-ph]]; V. K. Oikonomou, "The 
2d Gross-Neveu Model for Pseudovector Fermions and Tachyonic Mass Generation," 
arXiv:1109.6170 [hep-th]; P. Wang, H. Wu and H. Yang, "Superluminal neutrinos and 
domain walls," larXiv: 1109.69301 [hep-ph]; M. A. Anacleto, F. A. Brito and E. Pas- 
sos, "Supersonic Velocities in Noncommutative Acoustic Black Holes," arXiv:1109.6298 
[hep-th]. F. Giacosa, P. Kovacs and S. Lottini, "Could the OPERA setup send a bit 
of information faster than light?," larXiv: 11 10.36421 [hep-ph]. T. R. Morris, "Off-shell 
OPERA neutrinos," arXiv:1110.3266 [hep-ph]. J. Bramante, "Sterile Neutrino Pro- 
duction Through a Matter Effect Enhancement at Long Baselines," ar Xiv:1110. 4871 
[hep-ph]. 

[3] A. Mecozzi, M. Bellini, "Superluminal group velocity of neutrinos," 
|arXiv: 11 10. 12531 [hep-ph]]. 

[4] S. Tanimura, "Apparent Superluminal Muon-neutrino Velocity as a Manifestation of 
Weak Value," larXiv: 11 10. 17501 [hep-ph]. 



11 



[5] M. V. Berry, N. Brunner, S. Popescu and P. Shukla, "Can apparent superluminal 
neutrino speeds be explained as a quantum weak measurement?," J. Phys. A 44 (2011) 
492001 |arXiv:1110.2832l [hep-ph]]. 

[6] A. Nicolaidis, "Neutrino Shortcuts in Spacetime," arXiv:1109.6354 [hep-ph]. 

[7] Wang Yun-ping and Zhang Dian-lin, "Reshaping, path uncertainty, and superluminal 
traveling," Phys. Rev. A 52, 2597-2600 (1995). 

[8] M. M. Sanchez-Lopez, A. Sanchez-Meroho, J. Arias, J. A. Davis, and I. Moreno, "Ob- 
servation of superluminal and negative group velocities in a MachZehnder interferom- 
eter," Appl. Phys. Lett. 93, 074102 (2008). 

[9] Nicolas Brunner, Valerio Scarani, Mark Wegmiiller, Matthieu Legre, and Nicolas Gisin, 
"Direct Measurement of Superluminal Group Velocity and Signal Velocity in an Optical 
Fiber," Phys. Rev. Lett. 93, 203902 (2004). 

[10] K. Nakamura et al. (Particle Data Group), J. Phys. G37, 075021 (2010) and 2011 
partial update for the 2012 edition. 

[11] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo, A. M. Rotunno, "Evidence of 9i 3 > 
from global neutrino data analysis," 
|arXiv:1106.6028l [hep-ph]]. 

[12] A. E. Ball et al, "CNGS muon monitoring: Why do we need both muon detector 
stations at the start-up of the CNGS neutrino beam?," CERN SL Division tech- 
nical note SL-Note-2002-040 EA EMDS No. 36258 available on the CNGS website 



http://proj-cngs.web.cern.ch 



[13] N. Agafonova, A. Anokhina, S. Aoki, A. Ariga, T. Ariga, L. Arrabito, D. Autiero, 
A. Badertscher et al., "The Detection of neutrino interactions in the emulsion/lead 
target of the OPERA experiment," JINST 4 (2009) P06020. [arXiv:0903.2973 [hep- 
exl]. 



[14] http : / /cdsweb . cern. ch /record / 1 384486 



[15] N. Agafonova et al. [ OPERA Collaboration ], "Study of neutrino interactions with 
the electronic detectors of the OPERA experiment," New J. Phys. 13 (2011) 053051. 
|arXiv:1102.i882ll [hep-ex]] . 



12 



